Finite-temperature magnetism of transition metals: 
an LDA+DMFT approach. 
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We present an ab initio quantum theory of the finite temperature magnetism of iron and nickel. 
A recently developed technique which combines dynamical mean-field theory with realistic electronic 
structure methods, successfully describes the many-body features of the one electron spectra and 
the observed magnetic moments below and above the Curie temperature. 



The theory of itinerant electron ferromagnetism is one 
of the central problems in condensed matter physics (for 
reviews see There is a need for a first principles 

approach which is able to describe ground state and ther- 
modynamical quantities, as well as the one particle spec- 
tral properties of itinerant magnets. These quantities are 
currently being probed in spin polarized tunneling as well 
as spin polarized photoemission experiments with a view 
to possible applications as spin valves . 

Iron and nickel are the oldest and experimentally best 
studied prototypical systems, and serve as benchmarks 
for electronic structure methods. At very low temper- 
atures, a band like description of Fe and Ni has been 
very successful. The density functional theory (DFT) in 
the local density approximations (LDA) gives a quantita- 
tively accurate description of several ground state prop- 
erties of these materials such as the ordered magnetic 
moment and the spin wave stiffness [Q as calculated from 
the spin wave dispersion. 

While density functional theory, can in principle pro- 
vide a rigorous description of the finite temperature ther- 
modynamic properties, at present there is no accurate 
practical implemention available. As a result the finite 
temperature properties of magnetic materials are esti- 
mated following a simple suggestion jq], whereby con- 
strained DFT at T=0 is used to extract exchange con- 
stants for a classical Heisenberg model, which in turn 
is solved using approximation methods (e.g. RPA, mean 
field ) from classical statistical mechanics of spin systems 
The most recent implementation of this approach 
gives good values for the transition temperature of iron 
but not of nickel ||. While these localized spin models 
give, by construction, at high temperatures a Curie- Weiss 
like magnetic susceptibility, as observed experimentally 
in Fe and Ni, they encountered difficulties in predicting 
the correct values of the Curie constants [p~o|[ . 

It has been recognized for a long time, that to de- 
scribe the finite temperature aspects of itinerant elec- 
tron magnets, one needs a formalism that takes into ac- 
count the existence of local magnetic moments present 



above Tc ]TT| , p^ |. This one of the central problems in 
the physics of strongly correlated electron systems, which 
forces us to reconcile the dual character of the electron, 
which as a particle requires a real space description and 
as a wave requires a momentum space description in a 
unified framework. A very successful method satisfy- 
ing these requirements, the dynamical mean-field the- 
ory (DMFT) fl3[| , has been recently developed. This 
many-body scheme can be combined with standard LDA 
band structure calculations to include the effects of a 
realistic band structure QJl5|. Such "LDA+DMFT" 
approach has been successfully applied for the compu- 
tations of electronic structure and spin-wave spectrum 
of iron |LfflL7|. Nevertheless the most difficult and in- 
teresting finite-temperature effects were not considered 
previously and are the subject of this letter. 

Here we present realistic LDA+DMFT calculations 
of finite temperature magnetic properties of iron and 
nickel. A numerically exact Quantum Monte-Carlo 
(QMC) scheme is used for the solution the DMFT equa- 
tions. We find that a consistent first principles descrip- 
tion of the magnetic properties and of the one electron 
spectra of iron and nickel is possible, within an approach 
that makes only two essential approximations: the local- 
ity of the electron self energy and of the particle hole 
irreducible vertex . 

There have been numerous efforts to construct a many 
body theory of these materials, for example using the T- 
matrix theory jl8| , a self-consistent moment method [^9| , 
local 3-body equations J2(],|l| and the GW J22j approxi- 
mations. The LDA+DMFT approach goes beyond these 
works in the treatment of a realistic orbitally degenerate 
band structure, and in the treatment of the many body 
interactions. 

We start with the LDA Hamiltonian in the tight- 
binding orthogonal LMTO representation Hi^ff (k) [g3| , 
where m describes the orbital basis set containing 3<i-, 4s- 
and Ap- states, k runs over the Brillouin zone (BZ). The 
interactions are parameterized by a matrix of screened 
local Coulomb interactions U mm ' and and a matrix of 



f 



exchange constants J mm ', which are expressed in terms 
of two screened Hubbard parameters U and J, describ- 
ing the average Coulomb repulsion and the interatomic 
ferromagnetic exchange respectively. We use the values 
U = 2,3 (3, 0) eV for Fe (Ni) and the same value of the 
interatomic exchange, J = 0.9 eV for both Fc and Ni, a 
result of constrained LDA calculations [^4|,[l4 15 1. 

Dynamical mean-field theory maps the many-body 
system onto a multi-orbital quantum impurity, i.e. a 
set of local degrees of freedom in a bath described by 
the Weiss field function Q. The impurity action (here 
n ma = ctn<j c mcr and c(t) = [c TOCr (r)] is a vector of Grass- 
man variables) is given by: 
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It describes the spin, orbital, energy and temperature de- 
pendent interactions of particular magnetic 3ri-atom with 
the rest of the crystal and is used to compute the local 
Greens function matrix: 



G CT (r-r') = -| J D[c,c+]< 



e- s '"c(T)c+(T') (2) 



(Z is the partition function) and the impurity self energy 

The Weiss field function is required to obey the self 
consistency condition fl4| , [l5| ], which restores transla- 
tional invariance to the impurity model description: 

G a (cj n ) = J^K^n + m)i - n LDA (k) - K)]- 1 (3) 



fi is the chemical potential defined self-consistently 
through the total number of electrons, uj n — (2n + l)irT 
are the Matsubara frequencies for temperature T = (3~ 1 
(n = 0, ±1, ...) and a is the spin index. The local ma- 
trix X^ c is the sum of two terms, the impurity self en- 
ergy and a so-called "double counting " correction, Ed c 
which is meant to substract the average electron-electron 
interactions already included in the LDA Hamiltonian. 
For metallic systems we propose the general form of dc- 
correction: (w) = (u) - \Tr a Y. a (0). This is 
motivated by the fact that the static part of the corre- 
lation effects are already well described in the density 
functional theory. Only the d-part of the self-energy 
is presented in our calculations, therefore = for 
s- and p- states as well as for non-diagonal d — s,p 
contributions. In order to describe the finite tempera- 
ture ferromagnetism of transition metals we use the non 
spin-polarized LDA Hamiltonian H. LDA (k) and accumu- 
late all temperature-dependent spin-splittings in the self- 
energy matrix (ui n ). 

We use the impurity QMC scheme for the solution of 
the multiband DMFT equations (25). The Hirsch dis- 
crete Hubbard-Stratonovich transformations introduces 



(2M - 1)M auxiliary Ising fields S™ m , where M is the 
orbital degeneracy of the d-states and calculate G„ (r) by 
an exact integration of the fermion degrees of freedom in 
the functional integral (Eq.(|J)) |0|. In order to sample 
efficiently all the spin configurations in the multi-band 
QMC scheme, it is important to use "global" spin-flips: 

[Smm'\ ~ * [~S7nm' <y \ m addition to the local moves of the 
auxiliary fields. The number of QMC sweeps was of the 
order of 10 5 . A parallel version of the DMFT program 
was used to sample the 45 Ising fields for 3<i-orbitals. We 
used 256 k-points in the irreducible part of the BZ for the 
k integration. 10 to 20 DMFT iterations were sufficient 
to achieve convergence far from the Curie point. Due 
to the cubic symmetry of the bcc-Fe and fcc-Ni lattices 
the local Green function is diagonal in the basis of real 
spherical harmonics. The spectral functions for real fre- 
quencies were obtained from th e Q MC data by applying 
the maximum entropy method |26|. 
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FIG. 1. LDA+DMFT results for ferromagnetic iron 
(T = 0.8 Tq). The partial densities of d-states (full lines) is 
compared with the corresponding LSDA results at zero tem- 
perature (dashed lines) for the spin-up (red lines, arrow-up) 
and spin-down (blue lines, arrow-down) states. The insert 
shows the spin-spin autocorrelation function for T=1.2 Tc- 

Our results for the local spectral function for iron 
and nickel are shown in Figs.l and 2, respectively. The 
LDA+DMFT approach describes well all the qualitative 
features of the density of states (DOS), which is espe- 
cially non-trivial for nickel. Our QMC results reproduce 
well the three main correlation effects on the one par- 
ticle spectra below Tc ^7j: the presence of a famous 
6 eV satellite, the 30% narrowing of the occupied part 
of d-band and the 50% decrease of exchange splittings 
compared to the LDA results. Note that the satellite in 
Ni has substantially more spin-up contributions in agree- 
ment with photoemission spectra |^7j . Correlation effects 
in Fe are less pronounced than in Ni, due to its large spin- 
splitting and the characteristic bcc-structural dip in the 
density of states for spin-down states near Fermi level, 
which reduces the density of states for particle hole exci- 
tations. 
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FIG. 2. Same quantities as in Fig.l for ferromagnetic 
nickel (T = 0.9 Tc)- The insert shows the spin-spin autocor- 
relation function for T=1.8 Tc- 

The uniform spin susceptibility in the paramagnetic 
state: 



Xq=0 



dM 



(4) 



was extracted from the QMC simulations by measuring 
the induced magnetic moment in a small external mag- 
netic field. It includes the polarization of the impurity 
Weiss field by the external field E| . 

The dynamical mean field results accounts for the ob- 
served Curie -Weiss law which is observed experimentally 
in Fe and Ni. As the temperature increases above Tc, 
the atomic character of the system is partially restored 
resulting in an atomic like susceptibility with an effective 
moment: 



Xq=0 



fitf 



3(T-T C ) 



(5) 



The temperature dependence of the ordered magnetic 
moment below the Curie temperature and the inverse 
of the uniform susceptibility above the Curie point are 
plotted in Fig. 3 together with the corresponding experi- 
mental data for iron and nickel The LDA+DMFT 
calculations describes the magnetization curve and the 
slope of the high-temperature Curie- Weiss susceptibil- 
ity remarkably well. The calculated values of high- 
temperature magnetic moments extracted from Eq.(Q) 
are fi e ff = 3.09 (1.50)/iB for Fe (Ni), in good agreement 
with the experimental data fJ* e ff = 3.13 (1.62) /is for Fe 
(Ni) f§. 

We have estimated the values of the Curie tempera- 
tures of Fe and Ni from the disappearance of spin polar- 
ization in the self-consistent solution of DMFT problem 
and from the Curie- Weiss law in (Eq.(||)). Our estimates 
Tc = 1900 (700)K are in reasonable agreement with ex- 
perimental values of 1043 (631)K for Fe (Ni) respectively 
p8[ , considering the single site nature of the DMFT ap- 
proach. 



Within dynamical mean field theory one can also com- 
pute the local spin susceptibility defined by 



Xlo 
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dr (S(r)S(0)> 



(6) 



where g s — 2 is the gyromagnetic ratio and S = 
\ Y^ m ,a,<j> c lna a <™ ,c ma' is single-site spin operator and 
g = (o~ x , &y i c z ) are Pauli matrices. It differs from the 
<7 = susceptibility, Eq.(IJ), by the absence of spin polar- 
ization in the Weiss field of the impurity model. Eq.(||) 
cannot be probed directly in experiments but it is eas- 
ily computed in DMFT-QMC. Its behavior as function 
of temperature, gives a very intuitive picture of the de- 
gree of correlations in the system. In a weakly correlated 
system we expect Eq.(^|) to be nearly temperature in- 
dependent, while in a strongly correlated system we ex- 
pect a leading Curie- Weiss behavior at high temperatures 
Xiocai — ftfoc/fiT + const) where /i; oc is an effective local 
magnetic moment. In the Heisenberg model with spin S, 
/if oc = S(S+l)gg and for well-defined local magnetic mo- 
ments (e.g., for rare earth magnets) this quantity should 
be temperature independent. For the itinerant electron 
magnets \ii oc is temperature-dependent, due to a variety 
of competing many body effects such as Kondo screening, 
the induction of local magnetic moment by temperature 
|l2| and thermal fluctuations which disorders the mo- 
ments 1 29 . All these effects are included in the DMFT 
calculations. The r-dependence of the correlation func- 
tion (S (t) S(0)) results in the temperature dependence 
of ni oc and is displayed in the inserts on the Figs. 1,2. Iron 
can be considered as a magnet with very well-defined lo- 
cal moments above Tc (the T-dependence of the correla- 
tion function is relatively weak), whereas nickel is more 
itinerant electron magnet ( stronger r -dependence of the 
local spin- spin autocorrelation function ). 
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FIG. 3. Temperature dependence of ordered moment and 
the inverse ferromagnetic suscpetibility for Fe (open square) 
and Ni (open circle) compared with experimental results for 
Fe (square) and Ni (circle) (from Ref.[l,28]). 
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The comparison of the values of the local and the q = 
susceptibility gives a crude measure of the degree of short 
range order which is present above Tc- As expected, the 
moments extracted from the local susceptibility Eq.(|(]) 
are a bit smaller (2.8 fj,B for iron and 1.3 [Ib for nickel) 
than those extracted from the uniform mangetic suscep- 
tibility. This reflects the small degree of the short-range 
correlations which remain well above Tc |30[ . The high- 
temperature LDA+DMFT clearly show the presence of 
a local-moment above Tc- This moment, is correlated 
with the presence of high energy features (of the order of 
the Coulomb energies) in the photomeission. This is also 
true below Tc, where the spin dependence of the spec- 
tra is more pronounced for the satellite rigion in nickel 
than for that the quasiparticle bands near the Fermi level 
(Fig. 2). This can explain the apparent discrepancies be- 
tween different experimental determinations of the high- 
temperature magnetic splittings ||^J[lLp2| as being the 
result of probing different energy regions. The resonant 
photoemission experiments |3l| reflects the presence of 
local-moment polarization in the high-energy spectrum 
above Curie temperature in nickel, while the low-energy 
ARPES investigations |}2| results in non-magnetic bands 
near the Fermi level. This is exactly the DMFT view on 
the electronic structure of transition metals above Tc- 
Fluctuating moments and atomic-like configurations are 
large at short times, which results in correlation effects 
in the high-energy spectra such as spin-multiplct split- 
tings. The moment is reduced at longer time scales, cor- 
responding to a more band-like, less correlated electronic 
structure near the Fermi level. 

To conclude, we presented the first results of ab ini- 
tio LDA+DMFT calculations of finite-temperature mag- 
netic properties for Fe and Ni and showed that many 
body effects which incorporate the local atomic charac- 
ter of the electrons and which are ignored in the standard 
LDA based scheme are essential for a simultaneous de- 
scription of the magnetic properties and the one electron 
spectra of itinerant electron magnets. DMFT gives a 
satisfactory semiquantitative description of the physical 
propertis of Fe and Ni, far from the Curie point, indicat- 
ing that the critical fluctuations, which are not included 
in the DMFT approximation, do not play a crucial role, 
except for the immediate vicinity of the transition, and 
many aspects of the physics of this system can be under- 
stood within an approach which stresses local physics, ft 
would be interesting to extend this study to other itiner- 
ant magnetic systems with more atoms per unit cell, such 
as SrRuOa which is also well described by band theory 
at very low temperatures but has anomalous properties 
above its Curie point |33| . 
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